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Abstract 

Studies of linkage and association in various ethnic populations have revealed many predisposing genes of multiple 
neurotransmitter systems for alcohol use disorders (AUD). However, evidence often is contradictory regarding the 
contribution of most candidate genes to the susceptibility of AUD. We, therefore, performed a case-control study to 
investigate the possible associations of genes selected from multiple neurotransmitter systems with AUD in a 
homogeneous Tibetan community population in China. AUD cases (N = 281) with an alcohol use disorder identification 
test (AUDIT) score >10, as well as healthy controls (N = 277) with an AUDIT score <5, were recruited. All participants were 
genotyped for 366 single nucleotide polymorphisms (SNPs) of 34 genes selected from those involved in neurotransmitter 
systems. Association analyses were performed using PLINK version 1.07 software. Allelic analyses before adjustment for 
multiple tests showed that 15 polymorphisms within seven genes were associated with AUD (p<0.05). After adjustment for 
the number of SNPs genotyped within each gene, only the association of a single marker (rsl 0044881) in HTR4 remained 
statistically significant. Haplotype analysis for two SNPs in HTR4 (rsl 7777298 and rsl 0044881) showed that the haplotype 
AG was significantly associated with the protective effect for AUD. In conclusion, the present study discovered that the HTR4 
gene may play a marked role in the pathogenesis of AUD. In addition, this Tibetan population sample marginally replicated 
previous evidence regarding the associations of six genes in AUD. 
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Introduction 

Alcohol use disorders (AUD) — including alcohol dependence 
(AD) and alcohol abuse (AA) — are worldwide problems that 
induce complex clinical symptoms by damaging most organs in the 
human body, as well as the central and peripheral nervous systems 
[1,2]. The risk of AUD is explained by interactions among 
multiple environmental and genetic factors in which the ethanol 
metabolic system and most neurotransmitter systems are involved 
[3,4]. Although many study results support the effect of the genes 
encoding different alcohol and acetaldehyde dehydrogenase on 
risk of AUD [5,6], studies on the role of neurotransmitter genes 
have documented inconsistent results. For instance, although the 
associations of Taql-Al polymorphism of dopamine receptor D2 
(DRD2) gene polymorphisms with AUD have been repeatedly 
reported in many studies [7,8,9], they have not been reported in 
numerous others [10]. 

There are two subgroups of gamma-aminobutyric acid (GABA) 
receptors (the GABA A receptors and the GABA B receptors). They 
both belong to the GABAergic system, which is a chief inhibitory 
neurotransmitter in the central nervous system. For this system, 



genetic investigators have focused on GABA A receptors and have 
documented conflicting evidence of the contributing role of the 
GABAa receptors on chromosome 4p (GABRA1, GABRB1, and 
GABRG1) and 5q33-34 (GABRA1, GABRA6, GABRB2, and 
GABRG2) in the development of AUD [1 1,12,13]. The association 
of the susceptivity, severity, and treatment response of AUD with 
some polymorphisms of serotonin transporter (5-HTT) and 
serotonin (5-HT) receptor genes had also been reported by several 
initial studies, and as a result, similar research in different 
populations has increased in recent years, but results are not 
conclusive [8,14,15]. 

Although animal models and imaging studies have suggested 
that some 5-HT receptors (such as HTR1A, HTR1B, and HTR3) 
might also play roles in the pathogenesis of AUD [16,17,18], only 
a limited number of genetic association studies were performed, 
yielding mixed results [19,20]. Similarly, genes from other 
neurotransmitter systems — such as the neuropeptide Y [MPT] 
gene of neuropeptide Y system [21]; the neuronal nicotinic 
acetylcholine receptor beta 2 {CHRNB2) gene of cholinergic system 
[22]; period circadian clock 1 {PERI); and period circadian clock 2 
(PER2) of period family [23] — were also found to be associated 
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with AUD in some studies but failed to be replicated in other 
research [24]. In Chinese ethnic groups, relatively few studies have 
focused on the role of neurotransmitter genes in the pathogenesis 
of alcohol use disorders. Most of them studied populations in 
Taiwan, including Chinese Han and aboriginal minority groups 
[25-31], and a few studies recruited participants from Mainland 
Chinese populations [32,33,34]. These studies also documented 
inconclusive evidence for several genes in neurotransmitter 
systems, including DRD2, 5-HTT, and HTR1B [25,29,30,32]. 

Nonetheless, aforementioned association studies have suggested 
the implicated role of candidate genes from neurotransmitter 
systems, including dopaminergic, GABAergic, glutamatergic, 
opioid, cholinergic, and serotonergic systems. Some of these genes 
are located in or close to chromosomal regions (in chromosomes 1 , 
2, 4, 5, 7, 11, 13, and 17) indicated by linkage studies 
[35,36,37,38]. The inconsistent results for AUD studies were 
considered caused by multiple AUD vulnerable genes with small 
and diverse effect sizes in different populations [39]. Recently, 
several genome-wide association studies (GWAS) on AUD 
identified a number of modestly associated genes [40,41,42]. For 
example, Kendler et al. [43] found that the calcium-activated 
potassium channel subunit alpha- 1 (KCNMA1) gene was associated 
with alcohol dependence in European Americans. Nevertheless, 
the results of most of those studies proved inconclusive after the 
conservative statistical correction, which further confirmed that 
the effect size of each candidate gene is weak and may differ in 
different populations. 

To overcome the limitations of small-number-marker studies 
and the expensive cost of GWAS, it is necessary to choose a 
compromised alternative strategy. The new strategy is to select 
candidate genes from multiple neurotransmitter systems in a 
focused array based on the neurobiology of addiction and 
previously documented evidence. In the present study, we thus 
performed a case-control study to detect the associations with 
AUD of 34 genes selected from multiple neurotransmitter systems 
by using a set of samples collected from a homogeneous Tibetan 
population in China [44]. 

Methods 

Participants and Clinical Assessments 

We performed an epidemiological survey in 3171 individuals 
selected by the stratified cluster-random sampling method from 
four suburban counties of Lhasa City in Tibet [44] . The Chinese 
interview version of the Alcohol Use Disorders Identification Test 
(AUDIT) [45] — the reliability and validity of which had been 
proven good for application in this population [44] — was used to 
assess alcohol use disorders of respondents in the survey. Using 
DSM-IV diagnostic criteria according to SCID-P assessment as 
the gold criteria, we identified that the best cut-off point of 
AUDIT for diagnosing AUD in this population is 10, with both 
sensitivity and specificity >0.85, although the cut-off point for 
alcohol use disorders should be 6 if all cases need to be included 
(sensitivity was 1.00, but specificity was only 0.59) [44]. All 
respondents with an AUDIT score £ 1 0 (n — 4 1 9) were defined as 
AUD cases and invited to participate in molecular genetic 
research. Among those invited to participate in the study, 383 
respondents agreed to be involved and provided oral mucosa 
samples. Accordingly, 350 respondents matched by age, gender, 
and years of education, with AUDIT score £5, were appropriately 
recruited as healthy control subjects for the molecular genetic 
study [46,47]. After excluding poor-quality DNA samples which 
were not suitable for clip analysis, the present study included 281 
cases and 277 healthy controls. Statistically significant deviation 



with p<0.01 between case group and control group was shown in 
AUDIT score (mean±SD for case group: 17.84±6.12; control 
group: 0.96±1.54). There was no significant difference (p>0.05) 
between the two groups in terms of sex ratios (case group: 164 
males and 1 17 females; control group: 163 males and 1 14 females), 
age (case group: 47. 70± 12.19 years; control group: 47. 89± 13.73 
years), and years of education (case group: 2.42 ±2.64 years; 
control group: 2.1 1 ±2.74 years). None of the control subjects had 
a family history of AUD or a lifetime history of any other 
psychiatric disorders. This study was approved by the Ethics 
Committee of West China Hospital, Sichuan University. The 
written informed consent forms were obtained from all partici- 
pants. 

Selection and Genotyping of Candidate SNPs 

The criteria of candidate genes selection in this study were 
based on molecular pathway or previously reported evidence from 
linkage, association, or high-throughput analyses that were 
indicated to be related to AUD [22,43,48,49]. The 34 genes 
selected from several systems — including serotonin, dopamine, 
GABA, opioid, NMDA, and stress (HPA axis and neuropeptide Y) 
— were assessed in the study. Genomic regions containing 4 kb 
upstream and 4 kb downstream of each candidate gene were 
screened for tagger SNPs using the University of California Santa 
Cruz (UCSC) Genome Bioinformatics database (http: //genome. 
ucsc.edu) and the Broad Institute tagger tool (www.broadinstitute. 
org/mpg/tagger/server.html). After the linkage disequilibrium 
(LD) patterns were calculated (r2s0.8), selected tagger SNPs were 
assessed using the criterion of Illumina GoldenGate Custom 
Panels Design Score >0.6. The Illumina 384 GoldenGate Array 
Matrix (Illumina, Inc.) for 366 SNPs of 34 genes was designed and 
used in the present study (Table 1). 

Oral mucosa samples were obtained from all participants, and 
genomic DNA was extracted according to standard procedures. 
DNA concentrations were quantified by NanoDrop 8000 spec- 
trophotometer and normalized to 50 ng/uL for the chip. All 
samples (5 (J.L) were genotyped using the custom panel individu- 
ally, according to the manufacturer's protocol, on an Illumina 
BeadStation 500 G GoldenGate genotyping platform at Sichuan 
University. 

Data Analysis 

Allelic and haplotype association analyses were performed using 
PLINK version 1.07 software (http://pngu.mgh.harvard.edu/ 
~purcell/plink) for all genotyping data. Samples with genotyping 
rates and SNPs with minor allele frequency of ^0.05 were 
excluded; Hardy-Weinberg equilibrium (HWE) was tested 
(p>0.001). 

For haplotype association analysis, the LD was calculated 
between the pair of markers. Because each gene included in this 
study was selected using candidate-gene strategy rather than 
through random selection, statistical adjustments for multiple tests 
were performed for the number of SNPs genotyped within each 
gene [50,51]. 

Results 

Sixty-six of the 558 participants were excluded for low 
genotyping (MIND>0.05). After frequency and genotyping 
pruning, 336 SNPs remained, and the total genotyping rate in 
the remaining participants was 0.9728. Before adjusting for 
multiple tests, statistically significant deviations with p<0.05 
between the case group and control group were observed in 15 
polymorphisms within seven genes (Table 2). Among them, five 
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Table 1. Selection of candidate genes and SNPs. 







Systems 


Candidate genes 


Number of selected SNPs 


Dopamine system 


DRD2 


10 




DRD3 


5 




COMT 


12 




MAOA 


8 




MAOS 


6 




SLC6A3 


11 


Gamma-aminobutyric acid (GABA) system 


GA8RA1 


10 




GA8RA2 


10 




GA8RA4 


3 




GABRA5 


8 




GA8RA6 


4 




GA8RB2 


32 




GA8RB3 


24 




GABRG 1 


6 


GABRG2 18 




GABRG3 


33 


Serotonin system 


SLC6A4 


4 


HTR1A 1 




HTR1B 


6 




HTR2C 


4 




HTR3A 


6 




HTR3B 


7 




HTR3E 


2 




HTR4 


17 




HTR6 


7 




HTR7 


12 


Period family 


PERI 


3 




PER2 


16 


Neuropeptide Y 


NPY 


6 


Diacylglycerol kinases 


DGKZ 


2 


Ca2+-activated potassium channels 


KCNMA1 


53 


Adenosine receptor 


ADORA2A 


4 


Cholinergic system 


CHRNB2 


2 


Guanine nucleotide binding protein (G protein) 


GNAS 


14 



doi:1 0.1 371 /journal.pone.0080206.t001 



SNPs, including two (rsl7777298 and rsl0044881) with p<0.01, 
were in HTR4 gene, and one and two SNPs were in GABRA1 and 
GABRB2, respectively (all three genes were located in chromosome 
5q31-34). The other seven polymorphisms were in KCJVMA1, 
SLC6A4, COMT, and GABRG 1, located in chromosomes 10, 17, 
22, and 4, respectively. 

After adjustment for the number of SNPs within each gene, a 
single significant deviation (rs 10044881, adjusted p value 
= 2. 122 x 10 ) in HTR4 remained, which was observed in allelic 
frequencies of any examined single marker. 

We then nominated the two polymorphisms with p<0.01 in 
HTR4 (rsl7777298 and rs 10044881, D' = 0.69) to perform a 
haplotype association analysis. The result showed that the 
haplotype AG carried by 20.6% of AUD cases, versus 29.6% of 
the controls, was associated with the protective effect for AUD. 



The association remained statistically significant even after 
correction for multiple tests (Table 3). 

Discussion 

In this case-control study, we investigated markers selected from 
34 genes in neurotransmitter systems to evaluate the role of 
potential candidate genes in the pathogenesis of AUD in the 
Tibetan population. Some markers in HTR4 are detected to be 
associated with AUD, and it is especially noteworthy that the 
association remained statistically significant after adjustment for 
multiple tests — not only for the single locus, but also for the 
haplotype. This finding positively echoes the finding in a GWAS 
study in a Korean population, which suggested that four 
interesting genes — the cholinergic receptor/muscarinic 3 
[GHRM3] gene; the phosphodiesterase 4D /cyclic AMP-specific 
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Table 2. Allelic frequencies of the markers with unadjusted P<0.05. 





SNP 


Gene 


Chr 


A1/A2 


A1 Frequency 




P 


OR 


P' 










Case 


Control 








rs 10044881 


HTR4 


5 


G/A 


0.225 


0.318 


1.248X10" 3 * 


0.623 


2.122x10~ 2 ** 


rs4263535 


GABRA1 


5 


G/A 


0.444 


0.359 


6.992 x10~ 3 * 


1.429 


6.992 x10~ 2 


rsl 7777298 


HTR4 


5 


T/A 


0.333 


0.256 


9.034 x10~ 3 * 


1.452 


0.154 


rs655797 


KCNMA 1 


10 


G/A 


0.335 


0.258 


9.579 x10~ 3 * 


1.446 


0.508 


rsl 2782077 


KCNMA 1 


10 


A/G 


0.116 


0.174 


1.157x10~ 2 * 


0.625 


0.613 


rs7721 747 


HTR4 


5 


A/G 


0.263 


0.195 


1.209x10~ 2 * 


1.476 


0.206 


rs7712001 


HTR4 


5 


A/C 


0.337 


0.267 


1.912x10~ 2 * 


1.392 


0.325 


rs547212 


KCNMA 1 


10 


G/A 


0.121 


0.174 


2.021 x10~ 2 * 


0.651 


1.071 


rsl 0477387 


HTR4 


5 


G/C 


0.419 


0.348 


2.597 x10~ 2 * 


1.346 


0.441 


rs25528 


SLC6A4 


17 


C/A 


0.179 


0.237 


2.793 x10~ 2 * 


0.702 


0.112 


rs5993882 


COMT 


22 


C/A 


0.230 


0.174 


3.032 x10~ 2 * 


1.420 


0.303 


rsl 504501 


GA8RG1 


4 


G/A 


0.221 


0.281 


3.363 x10~ 2 * 


0.727 


0.202 


rsl 1080122 


SLC6A4 


17 


A/G 


0.172 


0.227 


3.398 x10~ 2 * 


0.708 


0.136 


rsl 7059393 


GABRB2 


5 


A/G 


0.185 


0.135 


3.511 x10~ 2 * 


1.454 


1.124 


rsl 7059409 


GABRB2 


5 


G/A 


0.187 


0.139 


4.450 x10~ 2 * 


1.425 


1.424 



Note: Chr, Chromosome; A, Allele; OR, odds ratio; 
p, unadjusted p values,* p<0.05; 
p', adjusted p values, ** p'<0.05. 
doi:1 0.1 371 /journal.pone.0080206.t002 



(PDE4D) gene; the neuronal periodic acid Schiff (PAS) domain 
protein 3 (NPAS3) gene; and the HTR4 gene — might be 
associated with amount of alcohol consumption from additive- 
model analyses with pooled data [52]. 

The human HTR4 gene was mapped to chromosome 5q3 1—33 
[53]. Its product is a glycosylated transmembrane protein that 
fulfills the role of modulating various neurotransmitters' release 
and enhancing synaptic transmission in the central nervous 
system. The gene thus had been proposed to be implicated in 
anxiety and cognitive function [54]. Several studies have shown 
that HTR4 polymorphisms could predispose to attention deficit 
hyperactivity disorder (ADHD) and bipolar disorder [55,56], while 
AUD was strongly associated with both ADHD and bipolar 
disorder [57,58,59,60]. Because the 5-HT4 receptor can modulate 
release of neurotransmitters (acetylcholine, dopamine, serotonin, 
and GABA) [6 1 ,62] , it is possible that the genetic polymorphism of 
HTR4 influences the interaction of neurotransmitters in the 
pathogenesis of AUD, which result in the association of AUD with 
polymorphisms of HTR4. 



In our study, in addition to the uniquely associated SNP, seven 
among the 1 4 SNPs marginally associated with AUD were located 
in the chromosome region 5q31-34, including four within HTR4 
and three within the GABAa receptor gene family (rs4263535 of 
GABRA1, rsl 7059393 and rsl 7059409 of GABRB2). Previous 
research has shown that four strong AUD candidate GABAa 
receptor genes(GABRAl, GABRA6, GABRB2, and GABRG2)- 
clustered together on locus 5q33-34 and are associated with AUD 
and its related traits [11,12,13]. It is thus possible that the 
association of markers in the HTR4 gene with AUD may reflect an 
indication on the chromosome region surrounding HTR4, rather 
than by HTR4 itself, due to linkage disequilibrium among markers 
on HTR4 and its surrounding region. 

Nevertheless, there is evidence against this possibility. For 
example, the distance from HTR4 to either GABRA1 or GABRB2 
was more than 10 Mb. Furthermore, our pair-wise LD analyses 
for rsl 0044881 (the unique associated SNP of HTR4) with 
rs4263535 (the marginally associated SNP of GABRA1), 
rsl 7059393, and rsl 7059409 (the two marginally associated SNPs 



Table 3. Estimated haplotype frequencies in the case-control subjects. 





Haplotype 


rsl 7777298 


rsl 0044881 


Frequency 




P 


P' 








Case 


Control 






1 


T 


G 


0.027 


0.019 


0.442 


>1 


2 


A 


G 


0.206 


0.296 


7.759 x10~ 4 * 


1.629x10~ 2 ** 


3 


T 


A 


0.307 


0.240 


0.014 


0.294 


4 


A 


A 


0.460 


0.445 


0.614 


>1 



Note: p, unadjusted p values,* p<0.05; 
p', adjusted p values, ** p'<0.05. 
doi:1 0.1 371 /journal.pone.0080206.t003 
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of GABRB2), respectively, showed very weak LD (D'^0.12) 
between them. 

Apart from the novel finding of the association between AUD 
and HTR4, this study using a Tibetan sample also marginally 
replicated previous evidence regarding the associations of six genes 
(GABRA1, GABRB2, KGNMA1, SLC6A4, COMT, and GABRG1) 
with AUD. In previous studies, all six of the genes were found to 
be associated with AUD and related traits, although revealed effect 
sizes all were small. For example, the GABA A receptor subunit 
genes on 5q33-34 have been found to be involved in psychological 
and behavioral disorders [63]. They demonstrated important roles 
in the acute and chronic effects of alcohol on the central nervous 
system [64] . In addition, various evidence has indicated that these 
four GABA A receptor genes may contribute to the pathogenesis of 
AUD in several populations, including that of East Asia 
[65,66,67]. For example, Park et al. [12] reported that polymor- 
phisms of the GABRA1 and GABRA6 receptor gene were associated 
with the development of alcoholism, and the polymorphisms of 
GABRA1 receptor were associated with the onset of alcoholism and 
alcohol withdrawal symptoms in the Korean population. Howev- 
er, several studies using the data from the Collaborative Study on 
the Genetics of Alcoholism (COGA) showed a weak relationship 
between the chromosome 5 cluster of GABA A receptor genes and 
AUD [13,68]. Nevertheless, more evidence has been presented on 
the association of chromosome 4 cluster of GABA A receptor genes 
with AUD and related traits from investigators [69] . Kendler et al. 
[43] performed a GWAS for the symptoms of alcohol dependence 
in the European American and African American participants and 
discovered the most significant intragenic SNP of KCNMA1 in the 
European American sample. Although results were inconsistent, 
many studies still suggested that the serotonin transporter gene (5- 
HTT) is an important candidate gene for AUD and related traits 
[70,14]. Extensive studies on a functional polymorphism (Vall58- 
Met) in the human catechol-O-mefhyltransferase [COMT] gene 
showed that the polymorphism (Vall58Met) was one of the 

References 

1. Baumberg B, Anderson P (2008) Trade anel health: how World Trade 
Organization (WTO) law affects alcohol and public health. Addiction 103: 
1952-1958. 

2. Testino G (2008) Alcoholic diseases in hepato-gastroenterology: a point of view. 
Hcpatogastroenterolog 55: 371—377. 

3. Schuckit MA. Smith TL (2006) An evaluation of the level of response to alcohol, 
externalizing symptoms, and depressive symptoms as predictors of alcoholism. 
J Stud Alcohol 67: 215-227. 

4. Timber-lake DS, Hopfcr CJ, Rhee SH, Friedman NP, Haberstick BC, et al. 
(2007) College attendance and its effect on drinking behaviors in a longitudinal 
study of adolescents. Alcohol Clin Exp Res 31: 1020-1030. 

5. Higuchi S, Matsushita S, Kashima H (2006) New findings on the genetic 
influences on alcohol nse and dependence. Curr Opin Psychiatry 19: 253-265. 

6. Zintzaras E, Stefanidis I, Santos M, Vidal E (2006) Do alcohol-metabolizing 
enzyme gene polymorphisms increase the risk of alcoholism and alcoholic liver 
disease? Hcpatology 43: 352-361. 

7. Blum K, Noble EP, Sheridan PJ, Montgomery A, Ritchie T, et al. (1990) Allelic 
association of humandopamine D2 receptor gene in alcoholism. J Am Med 
Assoc 263: 2055-2060. 

8. Foley PF, Loh EW, Innes DJ, Williams SM, Tannenberg AE, et al. (2004) 
Association studies of neurotransmitter gene polymorphisms in alcoholic 
Caucasians. Ann N Y Acad Sci 1025: 39-46. 

9. Noble EP, Zhang X, Ritchie TL, Sparkes RS (2000) Haplotypes at the DRD2 
locus and severe alcoholism. Am J Med Genet 96: 622—631. 

10. Heinz A, Sander T, Harms H, Finckh U, Kuhn S, et al.( 1996) Lack of allelic 
association of dopamine Dl and D2(TaqIA) receptor gene polymorphisms with 
reduced dopaminergic sensitivity to alcoholism. Alcohol Clin ExpRes 20:1109- 
1113. 

11. Enoch MA, Hodgkinson CA, Yuan Q, Albaugh B, Virkkuncn M, et al (2009) 
GABRG1 and GABRA2 as independent predictors for alcoholism in two 
populations. Neuropsyehopharmacology 34: 1245-1254. 

12. Park CS, Park SY, Lee CS, Sohn JW, Hahn GH, et al. (2006) Association 
between alcoholism and the genetic polymorphisms of the GABAA receptor 
genes on chromosome 5q33-34 in Korean population. J Korean Med Sci 21: 
533-538. 



significant markers of genetic predisposition to AUD [71,72]. 
COGA studies had linked alcohol dependence with the GABA A 
receptor gene family region that contains GABRG1 in the 
chromosome 4p [73,74]. Enoch et al. [11] found GABRG1 SNPs 
and haplotypes were significantly associated with AUD in both 
Plains Indians and Finnish Caucasians. In the present study, the 
marginal results of these genes suggest they may also play roles in 
the pathogenesis of AUD in the Tibetan population. Nevertheless, 
the statistical significance of the association between these genes 
and AUD was weakened after controlling for multiple tests. This 
may be due to the major limitation of the present study — i.e., the 
sample size is not large enough to detect those genes' small effect 
sizes after a conservative statistical correction. Therefore, our 
results need to be replicated in the larger and independent 
samples. 

In conclusion, the present study discovered that the HTR4 gene 
may play a marked role in the pathogenesis of AUD. In addition, it 
marginally replicated previous evidence regarding the associations 
of six genes (GABRA1, GABRB2, KCNMA1, SLC6A4, COMT, and 
GABRG1) with AUD in the Tibetan population. These findings 
enrich the knowledge regarding the genetic etiology of AUD. 

Acknowledgments 

We thank Jin Mei, Luo Dan, Sang Mu, Wang Mu, and Luosang Jiancai 
from Tibet University Medical College for their assistance with data 
collection. 

Author Contributions 

Conceived and designed the experiments: YX WJG TL XH XHL. 
Performed the experiments: YX YCW. Analyzed the data: YX QW XL. 
Contributed reagents/materials/ analysis tools: XHM WD QW. Wrote the 
paper: YX WJG TL XH. Collected the samples and information: WJG GL 
OL PZ BZ. 



13. Dick DM, Edcnbcrg HJ, Xuei X, Goate A, Hessclbrock V, et al. (2005) No 
association of the GABAA receptor genes on chromosome 5 with alcoholism in 
the collaborative study on the genetics of alcoholism sample. Am J Med 
Genet B Neuropsychiatr Genet 132B: 24-28. 

14. Fcinn R, Nellissery M, Kranzler HR (2005) Meta-analysis of the association of a 
functional serotonin transporter promoter polymorphism with alcohol depen- 
dence. Am J Med Genet B Neuropsychiatr Genet 133: 79-84. 

15. Dick DM, Plunkctt J, Hamlin D, Nurnberger J Jr, Kuperman S, et al. (2007) 
Association analyses of the serotonin transporter gene with lifetime depression 
and alcohol dependence in the Collaborative Study on the Genetics of 
Alcoholism (COGA) sample. Psychiatr Genet 17: 3.5-38. 

16. Martinez D, Slifstein M, Gil R, Hwang DR, Huang Y, et al. (2009) Positron 
emission tomography imaging of the serotonin transporterand 5-HT(lA) 
receptor in alcohol dependence. Biol Psychiatry 65: 175-180. 

1 7. Hu J, Henry S, GallezotJD, RopchanJ, NeumaierJF, et al. (2010) Serotonin IB 
Receptor Imaging in Alcohol Dependence. BIOL PSYCHIATRY 67: 800-803. 

18. Yoshimoto K, Yayama K, Sorimachi Y, Tani J, Ogata M, et al. (1996) 
Possibility of 5-HT3 receptor involvement in alcohol dependence: a micro- 
dialysis study of nucleus accumbens dopamine and serotonin release in rats with 
chronic alcohol consumption. Alcohol Clin Exp Res 20: 311A-319A. 

19. Polina ER, Contini V, Hutz MH, Bau CH (2009) The serotonin 2A receptor 
gene in alcohol dependence and tobacco smoking. Drug Alcohol Depend 101: 
128-131. 

20. Gorwood P, Lanfumey L, Hamon M (2004) Alcohol dependence and 
polymorphisms of serotonin-related genes. Med Sci (Paris) 20: 1132—1138. 

2 1 . Lappalainen J, Kranzler HR, Malison R, Price LH, Van Dyck C, et al. (2002) A 
functional neuropeptide Y Leu7Pro polymorphism associated with alcohol 
dependence in a large population sample from the United States. Arch Gen 
Psychiatry 59: 825-831. 

22. Ehrmgcr MA, Clegg HV, Collins AC, Corley RP, Crowley T, et al. (2007) 
Association of the neuronal nicotinic receptor bcta2 subunit gene {CIIRNB2) 
with subjective responses to alcohol and nicotine. Am J Med Genet B 
Neuropsychiatr Genet 144B: 596-604. 



PLOS ONE | www.plosone.org 



5 



November 2013 | Volume 8 | Issue 11 | e80206 



Genes Associated with Alcohol Use Disorders 



23. Dong L, Bilbao A, Laucht M, Henriksson R, Yakovlcva T, et al. (201 1) Effects of 
the circadian rhythm gene period 1 (perl) on psychosocial strcss-induccdalcohol 
drinking. Am J Psychiatry 168: 1090-1098. 

24. Kochnke MD, Schick S, Lutz U, Willccke M, Koehnke AM, et al. (2002) 
Severity of Alcohol Withdrawal Symptoms and the T1128C Polymorphism of 
the Neuropeptide Y Gene. J Neural Transm 109: 1423-1429. 

25. Lee JF, Lu RB, Ko HC, Chang FM, Yin SJ, Pakstis AJ, et al. (1999) No 
association between DRD2 locus and alcoholism after controlling the ADH and 
ALDH genotypes in Chinese Han population. Alcohol Clin Exp Res 23:592- 
599. 

26. Chen WJ, Chen CH, Huang J, Hsu YP, Scow SV, Chen CC, et al. (2001) 
Genetic polymorphisms of the promoter region of dopamine D2 receptor and 
dopamine transporter genes and alcoholism among four aboriginal groups and 
HanChincse in Taiwan. Psychiatr Genet 11:187-195. 

27. Loh cl-W, Farm CS, Chang YT, Chang CJ, Cheng AT (2004) Endogenous 
opioid receptor genes and alcohol dependence among Taiwanese Han. Alcohol 
Clin Exp Res:15-19. 

28. Loh cl-W, Lane HY, Chen CH, Chang PS, Ku LW, et al. (2006) Glutamate 
decarboxylase genes and alcoholism in Han Taiwanese men. Alcohol Clin Exp 
Res 30:1817-1823. 

29. Huang SY, Lin WW, Wan FJ, Chang AJ, Ko HC, et al. (2007) Monoamine 
oxidase-A polymorphisms might modify the association between the dopamine 
D2 receptor gene and alcohol dependence. J Psychiatry Neurosci 32: 185-192. 

30. Lcc SY, Lin WW, Huang SY, Kuo PH, Wang CL, et al.(2009) The relationship 
between serotonin receptor IB polymorphisms A- 16 IT and alcohol dependence. 
Alcohol Clin Exp Res 33:1589-1595. 

31. Lee SY, Hahn CY, Lee JF, Huang SY, Chen SL, et al. (2010) MAOA interacts 
with the ALDH2 gene in anxiety-depression alcohol dependence. Alcohol Clin 
Exp Res 34:1212-1218. 

32. Cao JX, HuJ, Yc XM, Xia Y, Haile CA, Kosten TR, et al. (201 1) Association 
between the 5-HTR1B gene polymorphisms and alcohol dependence in a Han 
Chinese population. Brain Res 1376:1—9. 

33. Wang XJ, Zhong SR, Bao JJ, Don SJ, Wu WY, et al. (2011) Association of 
polymorphism in the serotonin transporter gene promote with the susceptibility 
to alcohol dependence in Yunnan Han Population. Yi Chuan 33:48—53. 

34. Huang C, Jin AH, Li CH, Han HM, Piao XX (2012) The relation between 
OPRM1 gene A118G polymorphism and alcoholic liver disease, alcohol- 
dependent in Korean-Chinese and Han nationality men of Yanbian area. 
Journal of Medical Science Yanbian University 35:199—205. 

35. Agrawal A, Hinrichs AL, Dunn G, Bertelsen S, Dick DM, et al. (2008) Linkage 
scan for quantitative traits identifies new regions of interest for substance 
dependence in the Collaborative Study on the Genetics of Alcoholism (COGA) 
sample. Drug Alcohol Depend 93: 12-20. 

36. Gclernter J, Kranzlcr HR, Panhuysen C, Weiss RD, Brady K, et al. (2009) 
Dense genomewide linkage scan for alcohol dependence in African Americans: 
significant linkage on chromosome 10. Biol Psychiatiy 65: 111—115. 

37. Hill SY, Shen S, Zezza N, Hoffman EK, Perlin M, et al. (2004) A genome wide 
search for alcoholism susceptibility genes. Am J Med Genet B Neuropsychiatr 
Genet 128B: 102-113. 

38. Kuo PH, Nealc MC, Riley BP, Webb BT, SuUivan PF, ct al. (2006) 
Identification of susceptibility loci for alcohol-related traits in the Irish Affected 
Sib Pair Study of Alcohol Dependence. Alcohol Clin Exp Res 30: 1807-1816. 

39. Nurnbergcr JI Jr, Bicrut LJ (2007) Seeking the connections: alcoholism and our 
genes. Sci Am 296: 46-53. 

40. Bierut LJ, Agrawal A, Bucholz KK, Dohcny KF, Pugh E, et al. (2010) A 
genomc-wide association study of alcohol dependence. Proc Nad Acad Sci USA 
107: 5082-5087. 

41. Edcnberg HJ, Koller DL, Xuci X, Wctherill L, McClintick JN, ct al. (2010) 
Genome-wide association study of alcohol dependence implicates a region on 
chromosome 11. Alcohol Clin Exp Res 34: 840-852. 

42. Treutlcin J, Cichon S, Ridinger M, Wodarz N, Soyka M, et al. (2009) Genome- 
wide association study of alcohol dependence. Arch Gen Psychiatry 66: 773- 
784. 

43. Kcndler KS, Kalsi G, Holmans PA, Sanders AR, Aggen SH, ct al. (2011) 
Genomewide association analysis of symptoms of alcohol dependence in the 
molecular genetics of schizophrenia (MGS2) control sample. Alcohol Clin Exp 
Res. 35: 963-975. 

44. Guo W, Lanzi G, Luobu O, Ma X, Zhcn P, et al. (2008) An epidemiological 
survey of alcohol use disorders in a Tibetan population. Psychiatry Res 159: 56— 

66. 

45. Babor TF, John C, Saunders JB, Maristela G (2001) AUDIT: The Alcohol Use 
Disorders Identification Test: Guidelines for Use in Primary Health Care. 2nd 
edition. Geneva, Switzerland: World Health Organization/MSD/MSB. 

46. Guo W, Wang Q, Lanzi G, Luobu O, Ma X, Wang Y, et al. (2010) Interaction 
among genes influencing cthanol metabolism and sex is association with alcohol 
use disorders in a Tibet population. Am J Med Genet B Neuropsychiatr Genet 
153B:561-9. 

47. Guo WJ, Wang Q, Lanzi G, Luobu O, Ma XH, Wang YC, et al. (201 1) Gender- 
specific interactions between alcohol metabolism genes and severity of 
quantitative alcohol-rclatcd-traitsin a Tibetan population. Neurosci Lett 
495:22-5. 



48. Kohnkc MD (2008) Approach to the genetics of alcoholism: a review based on 
pathophysiology. Biochcm Pharmacol 75:160-177. 

49. Kimura M, Higuchi S (2011) Genetics of alcohol dependence. Psychiatry Clin 
Neurosci 65:213-225. 

50. Lan Q, Morton LM, Armstrong B, Hartgc P, Menashc I, et al. (2009) Genetic 
variation in caspase genes and risk of non-Hodgkin lymphoma: a pooled analysis 
of 3 population-based case-control studies. Blood 114:264-267. 

51. Perlis RH, Fijal B, Dharia S, Houston JP (2013) Pharmacogcnctic investigation 
of response to duloxetinc treatment in generalized anxiety disorder. 
PharmacogenomicsJ 13:280-285. 

52. Baik I, Cho NH, Kim SH, Han BG, Shin C (201 1) Genome-wide association 
studies identify genetic loci related to alcohol consumption in Korean men. 
Am J Clin Nutr 93: 809-16. 

53. Claeysen S, Faye P, Sebbcn M, Lcmairc S, BockacrtJ, ct al. (1997) Cloning and 
expression of human 5-HT4S receptors. Effect of receptor density on their 
coupling to adcnylyl cyclase. NcuroRcport 8: 3189-3196. 

54. King MV, Marsden CA, Fone KCF (2008) A role for the 5-HT(lA), 5-HT4 and 
5-HT6 receptors in learning and memory. Trends Pharmacol Sci 29: 482-492. 

55. Li J, Wang Y, Zhou R, Wang B, Zhang H, et al. (2006) Association of attention- 
deficit /hyper activity disorder with serotonin 4 receptor gene polymorphisms in 
Han Chinese subjects. Neurosci. Lett 401:6—9. 

56. Ohtsuki T, Ishiguro H, Detera-Wadleigh SD, Toyota T, Shimizu H, et al. (2002) 
Association between serotonin 4 receptor gene polymorphisms and bipolar 
disorder in Japanese case-control samples and the NIMH Genetics Initiative 
Bipolar Pedigrees. Mol Psychiatry 7: 954-961. 

57. August GJ, Winters KC, Realmuto GM, Fahnhorst T, Botzct A, ct al. (2006) 
Prospective study of adolescent drug use among community samples of ADHD 
and non-ADHD participants. J Am Acad Child Adolesc Psychiatry 45: 824-832. 

58. Molina BS, Pelham WE, Gnagy EM, Thompson AL, Marshal MP (2007) 
Attention-deficit /hyperactivity disorder risk for heavy drinking and alcohol use 
disorder is age specific. Alcohol Clin Exp Res 31: 643—654. 

59. Jacob T, Waterman B, Heath A, True W, Bucholz KK, et al. (2003) Genetic and 
environmental effects on offspring alcoholism: new insights using an offspring-of- 
twins design. Archives of General Psychiatry 60:1265—1272- 

60. Haydcn EP, Nurnbergcr JI Jr (2006) Molecular genetics of bipolar disorder. 
Genes, Brain, and Behavior 5: 85-95. 

61. Bianchi C, Rodi D, Marino S, Beam L, Siniscalchi A (2002) Dual effects of 5- 
HT4 receptor activation on GABA release from guinea pig hippocampal slices. 
Neuroreport 13: 2177-2180. 

62. Fink KB, Gothert M (2007) 5-HT receptor regulation of neurotransmitter 
release. Pharmacol. Rev 59: 360—417. 

63. Kerner B, Brugman DL, Frcimcr NB (2007) Evidence of linkage to psychosis on 
chromosome 5q33-34 in pedigrees ascertained for bipolar disorder. Am J Med 
Genet B Neuropsychiatr Genet 144B: 74-78. 

64. Grobin AC, Matthews DB, Devaud LL, Morrow AL (1998) The role of GABAA 
receptors in the acute and chronic effects of cdianol. Psychopharmacology 139: 
2-19. 

65. Loh EW, Smith I, Murray R, McLaughlin M, McNulty S, et al. (1999) 
Association between variants at the GABAAbcta2, GABAAalpha6 and 
GABAAgamma2 gene cluster and alcohol dependence in a Scottish population. 
Mol Psychiatry 4: 539-544. 

66. Loh EW, Higuchi S, Matsushita S, Murray R, Chen CK, ct al. (2000) 
Association analysis of the GABA(A) receptor subunit genes cluster on 5q33-34 
and alcohol dependence in ajapanese population. Mol Psychiatry 5: 301-307. 

67. Radel M, Vallejo RL, Iwata N, Aragon R, LongJC, et al. (2005) Haplotype- 
bascd localization of an alcohol dependence gene to the 5q34 {gamma}- 
aminobutyric acid type A gene cluster. Arch Gen Psychiatiy 62: 47—55. 

68. Song J, Koller DL, Foroud T, Carr K, Zhao J, ct al. (2003) Association of 
GABA(A) receptors and alcohol dependence and the effects of genetic 
imprinting. Am J Med Genet B Neuropsychiatr Genet 117B: 39—45. 

69. Agrawal A, Edenberg HJ, Foroud T, Bicrut LJ, Dunne G, ct al. (2006) 
Association of GABRA2 with drug dependence in the collaborative study of the 
genetics of alcoholism sample. Behav Genet 36: 640—650. 

70. Konishi T, Calvillo M, Lcng A, Lin K, Wan Y (2004) Polymorphisms of the 
dopamine D2 receptor, serotonin transporter, and GABAA receptor b3 subunit 
genes and alcoholism in Mexican- Americans. Alcohol 32: 45—52. 

71. Kibitov AO, Voskobocva EIu, Brodianskii VM, Chuprova NA, Smirnova EV 
(2010) Association study of the Vall58Mct polymorphism of the catcchol-O- 
mc thy ltransfe rase gene and alcoholism and heroin dependence: the role of a 
family history. Zh Nevrol Psikhiatr Im S S Korsakova 1 10: 84-88. 

72. Kohnke MD, Wiatr G, Kolb W, Kdhnke AM, Schick S, et al. (2003) Plasma 
homovanillic acid: a significant association with alcoholism is independent of a 
functional polymorphism of the human catcchol-O-mcfhyl transferase gene. 
Neuropsychopharmacology 28:1 004- 1010. 

73. Reich T, Edenberg HJ, Goatc A, Williams JT, Rice JP, et al. (1998) Genome- 
wide search for genes affecting the risk for alcohol dependence. Am J Med Genet 
81: 207-215. 

74. Long JC, Knowler WC, Hanson RL, Robin RW, Urbanek M, et al. (1998) 
Evidence for genetic linkage to alcohol dependence on chromosomes 4 and 1 1 
from an autosomc-widc scan in an American Indian population. Am J Med 
Genet 81: 216-221. 



PLOS ONE | www.plosone.org 



6 



November 2013 | Volume 8 | Issue 11 | e80206 



